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It is believed that one can extract more accurate information of the pion distribution amplitude 
from the pion-photon transition form factor (TFF) due to the single pion in this process. However 
the BABAR and Belle data of the pion-photon TFF have a big difference for Q 2 G [15,40] GeV 2 , 
and at present, the pion DA can not be definitely determined from the pion-photon TFF. it is 
crucial to find the right pion DA behavior and to determine which data is more reliable. In this 
letter, we perform a combined analysis of the most existing data of the processes involving pion 
by using a general model for the pion wavefunction/DA. Such a DA model can mimic all the 
existed pion DA behaviors, whose parameters can be fixed by the constraints from the processes 
7T° — > 77, 7r — > fiu, B — > -kIv, and etc. Especially, we examine the B — > it transition form 
factors that provides another constraint to the parameter B in our DA model, which results 
in B £ [0.00,0.29]. This inversely shows that the predicted curve for the pion-photon TFF is 
between the BABAR and Belle data in the region Q 2 G [15,40] GeV 2 . It will be tested by coming 
more accurate data at large Q 2 region, and the definite behavior of pion DA can be concluded finally. 



PACS numbers: 12.38.-t, 12.38.Bx, 14.40.Aq 



The recent data by the Belle Collaboration [l| is dra- 
matically different from those reported by the BABAR 
Collaboration 0], i.e. instead of a pronounced growth 
of the TFF at high Q 2 region, observed by BABAR, the 
Belle data are compatible with the well-known asymp- 
totic prediction Q 2 F 7T1 (Q 2 )\q2^ oc — > y/2f v , where the 
pion decay constant f w = 130.41 ± 0.03 ± 0.20 MeV Q. 
One thinks that it would be a big challenge to QCD if the 
the BABAR data are confirmed. One should ask a ques- 
tion what is the realistic light-cone wave function of the 
pion. Theoretically the pion-photon TFF depends heav- 
ily on the pion DA and sensitive to it in the higher energy 
region. Such large discrepancy between the Belle and 
the BABAR data at the high Q 2 region, shows that the 
pion-photon TFF can't discriminate the different pion 
DA models. On other hand, All the processes involving 
pion can provide more strong constraints on the pion DA 
for determining it. It would be helpful to have a gen- 
eral pion DA model that can mimic all the DA behaviors 



suggested in the literature. 

Recently, we have raised a pion wave function (WF) 0, 
jfj following the idea of Refs.[(l-[8j], in which the BHL pre- 
scription Q and the Melosh rotation [l(| are applied. 
More explicitly, the full form of the pion WF is written 

as 



*0z,k ± ) = J2 X AlA2 (a-,k±)* fl (x,k ± ), 



(1) 



where \ Xl X2 (x, kj_) stands for the spin-space WF and the 
spatial WF 
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The pion DA can be derived by integrating out the trans- 
verse momentum dependence of the WF, and we obtain 
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where Erf (a;) = J* e * dt, /io stands for the factor- distribution broadness of the WF, which can be expanded 
ization scale. The^pion DA at any scale /i 2 can be de- in the Gegenbauer polynomials, and by keeping its first 
rived from the initial DA /ig) through QCD evolu- two terms, we have (^(x) = 1 + B x C 2 3 ^ 2 (2a; — 1) . 

tion [ll[. The function <pn(x) dominates the longitudinal 

In the above model, there are four parameters A, (3, 
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m q and B. All the processes involving pion can provide 
strong constraints on the pion DA. Then these parame- 
ters can be fixed by a combination of various constraints, 

• The process 7r — > [iv provides the WF normaliza- 
tion condition 



dx 
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• The process 7r° 
WF 01, 



77 provides a sum rule for the 

(5) 



[ dar*0,kj_ = 0) = V6/f n 
Jo 



• The requirement of a reasonable value for the light 
valence quark mass m q provides another constraint, 
whose value should be around 0.30 GeV. 

• The parameter B can be regarded as an effective 
parameter that absorbs the contributions from all 
higher Gegenbauer polynomials, which dominantly 
determines the broadness of the WF. Its value can 
be fixed by high energy processes involving pion. 



TABLE I: Typical pion WF parameters under the condition of 
m q — 0.30 GeV, and its corresponding Gegenbauer moments. 
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0.00 


25.01 


0.586 


0.027 


-0.027 


0.30 


20.22 


0.668 


0.362 


-0.018 


0.60 


16.59 


0.745 


0.679 


0.020 



It is obvious that one can determine the right pion 
DA behavior by comparing the theoretical predictions 
derived by using different pion DA with the experimen- 
tal data. The pion DA as shown by Eq.Q can be con- 
veniently adopted for such purpose. We compare this 
model with other existed models in the literature and 
find it can mimic all the existed pion DA behaviors, e.g. 
as shown in Table HI keeping m q = 0.30 GeV and varying 
the parameter B within the region of [0.00, 0.60], the pion 
DA ([3]) shall vary from asymptotic-like to CZ-like. A de- 
tailed discussion on this point can be found in Refs.[l,Q. 
Then, to determine the right DA behavior is equivalent 
to find a proper B that can well explain the data. 

The process 77* — > ir° provides the simplest exam- 
ple of the perturbative application to exclusive processes 
since this process relates two photons with a single pion. 
It is usually regarded as a suitable platform to determine 
the pion DA. However, the large discrepancy between the 
Belle and the BABAR data at the high Q 2 region doesn't 
help to discriminate the different pion DA models or to fix 
the parameter B. In fact, by varying B within the region 
of [0.00, 0.60], it has been shown that the estimated form 
factors at large Q 2 region are quite different and they are 



in between the BABAR and Belle data for various value 
of B 0,I[E1. 

In the literature, the B — > 7r transition form factor 
fB^nv j g usua jiy adopted to extract the CKM matrix el- 
ement \V u b\- A consistent analysis of fS ^ n within the 
whole physical region has been in Ref . [13[ . In the frame- 
work of QCD light-cone sum rules, the expression for 
jB->7r h eav y depends on the correlation function: dif- 
ferent choice of the currents in the correlation function 
shall result in different expressions, in which, the pionic 
different twist structures provide different contributions. 
More explicitly, as a direct estimation up to twist-4 level, 
the ff _i-7r contains the full pion twist-2, 3 and 4 contribu- 
tions by taking the conventional correlator On 
the other hand, by taking the correlator with a chiral cur- 
rent (l9l - [23l ] , the most uncertain twist-3 contributions can 
be eliminated; while by taking the correlator suggested 
by Refs.[13, [25|, only the twist-3 contributions are kept 
up to twist-4 level. 

In Refs. [22I l23l| . it has been pointed out that the pionic 
leading twist contribution is dominant, which provides a 
greater than 95% contribution to the form factors. This 
provides us a good change to study the properties of the 
leading twist pion DA, that is, by using the experimen- 
tal values for \V u b\ and f^ 7 ' {G)\V u b\, one can expect to 
obtain a strong constraint for the parameter B. 

On the experimental side, from the measured shapes 
of the form factor for B — > ttIv, the CKM matrix element 
I Kb I multiplied by ff^ n (0) is H^: 



+ 



T (0)|K & | = (9.4 ±0.3 ±0.3) 



x 10 
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The CKM matrix element |V„ 6 | is (3.23 ± 0.31) x 10 
which is derived from a simultaneous fit to the experi- 
mental partial rates and lattice points on the exclusive 
process B — } ttW versus q 2 27 1. As a combination, we 

i.e. 

(7) 



obtain an moderate experimental limit for /^"^(O), 
/*-*(0) = 0.291™. 



It is noted that this upper limit f+~* n (0) = 0.301 agrees 
with the theoretical limit derived in most of the litera- 
ture, either under the pQCD approach (cf. f+~**(0) = 
0.265 ± 0.032 [H]) or the QCD light-cone sum rules 
approach (cf. /f^O) = 0.258 ± 0.031 0) or the 
extrapolated lattice QCD approach (cf. ff^^iO) — 
0.27±0.11 [H). While the lowest limit 0.282 for ff^(0) 
is somewhat too strict, however a better determination 
on this value is less important here, since it only de- 
termines the lower limit for the parameter B, which al- 
ready reaches its lowest value B min = 0.00 for the present 
choice. 

On the theoretical side, we adopt the light-cone sum 
rules derived in Ref.[23| for our analysis. During the 
calculation, we take the present DA model ([3]) to do the 
calculation; while the other input parameters are chosen 
to be the same as those of Ref. [23| . We present the value 
°f /+^ T (0) versus the parameter B in Fig.([T]), where 
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FIG. 1: The value of f+~* n (0) versus the parameter B. The 
solid, the dashed and the dash-dot lines stand for central, 
the upper and the lower values for /+ _>7r (0) that are derived 
from the QCD light-cone sum rules [23|. The shaded band 
represents the value derived from a combination of the data 
on f^(0)\V ub \ and \V ub \ [MEl- 



As a summary, we have performed a combined anal- 
ysis of the most existing data of the processes involving 
pion by using a flexible model for the pion DA. Such a 
DA model can mimic all the existed pion DA behaviors, 
whose parameters can be fixed by the constraints from 
the processes 7r° — > 77, ir — > [iv, B — > nW, and etc. Es- 
pecially, the experimental data on the B — ¥ it transition 
form factor provides useful information for determining 
the pion DA. 

Using the experiments on the values of /+ at q 2 = 
GeV 2 , we show that the reasonable range of the effec- 
tive parameter B of pion DA is [0.00,0.29]. It is noted 
that such region is consistent with Brodsky and Tera- 
mond's holographic DA model whose second Gegenbaulcr 
moment is about 0.145 [3(j. Therefore all of parameters 
in our DA model are fixed to a certain confidence level. 
A better understanding of other pionic processes, such as 
D — > nlD, shall provide further constraints on the wave 
function parameters. 



the theoretical errors are the squared average of those 
dominant errors caused by taking the MS 6-quark mass 
m b (fh b ) = 4.164 ± 0.025 GeV [29|, the B meson decay 
constant up to ncxt-to-leading order fg = 214^ [18[ 

and the effective threshold sfi = (34 ± 0.5) GeV 2 . We 
obtain the range of the parameter B: 

B = [0.00, 0.29]. (8) 
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FIG. 2: 
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Q 2 F 7T ~ / (Q 2 ) by taking m q = 0.30 GeV and B G 



[0.00, 0.29], where as a comparison, the Belle jj, BABAR @] 
CELLO |3j] and CLEO [12l data are presented. 



The pion-photon TFF can be calculated exactly by 
using the pion DA light-cone wave function ([3]). In the 
small Q 2 region, Q 2 < 15 GeV 2 , the pion-photon TFF 
for B ~ [0.00, 0.29] can explain the Belle BABAR 0, 
CELLO [H| and CLEO [32[ experimental data. While 
for the large Q 2 region, as can be seen from Fig.©, the 
pion-photon TFF is sensitive to the parameter B. The 
result with B <G [0.00, 0.29] is consistent with previous ob- 
servation that a moderate pseudo-scalar (jr, 77 or rj') DA 
B ~ 0.1 — 0.3 that is close to asymptotic-like can roughly 
explain the TFFsQ 2 F 7r7 , Q 2 F, n and Q 2 F v i~ 1 data simul- 
taneously [1, [33|,|34[. It has shown from Fig.© that the 
curve with the parameter B — 0.00 is more close to the 
Belle data and the curve with the parameter B = 0.29 is 
between the BABAR and Belle experimental data. Thus, 
it will be tested by coming more accurate data at large 
Q 2 region, and the definite behavior of pion DA can be 
concluded finally. 
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